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We report a clear evidence of atomic fractals in the nonlinear motion of a two-level atom in a 
standing-wave microcavity. Fractal-like structures, typical for chaotic scattering, are numerically 
found in the dependencies of outgoing positions and momenta of scattered atoms on their ingoing 
values and in the dependence of exit times of cold atoms on their initial momenta in the generic 
semiclassical models of cavity QED (1) with atoms in a far-detuned amplitude (phase)-modulated 
standing wave and (2) with coupled atomic external and internal degrees of freedom. Tiny interplay 
between all the degrees of freedom in the second model is responsible for trapping atoms even in 
a very short microcavity. It may lead simultaneously, at least, to two kinds of atomic fractals, 
a countable fractal (a set of initial momenta generating separatrix-like atomic trajectories) and a 
seemingly uncountable fractal with a set of momenta generating infinite walkings of atoms inside 
the cavity. 

PACS numbers: 42.50.Vk, 05.45.Df 



1. Cavity quantum electrodynamics (QED) with cold 
atoms is a rapidly growing field of atomic physics and 
quantum optics studying the atom-photon interaction in 
cavities |l]] . Recent experiments |^, || succeeded in explo- 
ration coupled external atomic center-of-mass, internal 
atomic, and field dynamics under condition of strong- 
coupling between a single cold atom and a single-mode 
field in a high-finesse optical microcavity. Methods to 
monitoring the atomic motion in real time have been re- 
alized experimentally. They open a new way to study 
the very foundations of the matter-field interaction and 
use particles trapped within a cavity mode for quantum 
communications, to monitoring of chemical and biolog- 
ical processes at the single-molecule scale and for other 
purposes. 

The basic model Hamiltonian of the interaction of a 
two-level atom with a single-mode standing-wave field in 
an ideal cavity is given by 

—f^of(t)(a/a^+aa + )cos(kfx), (1) 

where x and p are the atomic position and momentum op- 
erators, a the atomic Pauli operators, a and a' the field- 
mode operators. We incorporate a function of time f(t) 
in (0) to describe a possible modulation of the standing- 
wave amplitude in the framework of the basic Hamilto- 
nian. The strongly coupled atom-field system (Q) is a 
highly nonlinear one and is known to possess a rich va- 
riety of qualitatively different dynamics including chaos. 
The full Hamiltonian (|l|) for an atom, to be placed in a 
far-detuned modulated standing-wave, can be reduced to 
an effective non-autonomous Hamiltonian of a "sructure- 
less" atom with one and half degree of freedom which is a 
paradigma for quantum chaos in atomic optics jj, H 0, Q ■ 
It has been shown in || ^ that semiclassical chaos in the 
sense of extremal sensitivity to small changes in initial 
conditions is possible without any modulation, i.e. with 



an autonomous dynamical system with three degrees of 
freedom, the field, external, and internal atomic ones. 
In this Letter we report with both the models a clear 
evidence of a new property of the atom- field dynamics, 
atomic fractals in cavity QED. In the strong-coupling 
limit, one neglects dissipation in all the degrees of free- 
dom that may be justified by available microcavities with 
very large quality factors Q > 10 6 || or/and operating 
with far-detuned atoms. 

2. Scattering of atoms by a far-detuned mod- 
ulated standing wave. Following a scheme of experi- 
ments on scattering of atoms by light Jl(| , suppose that 
a monokinetic atomic beam propagates at an angle to 
the direction y which is perpendicular to the standing- 
wave axis x. The field is assumed to be uniform in the 
y-dircction. In a reference frame moving with a constant 
velocity in the y-direction there remains only the trans- 
verse atomic motion which can be considered classically 
if the atomic momentum is much larger than the photon 
momentum. So the atom is treated as a point parti- 
cle subject to a position-dependent force and is thus de- 
flected. One measures atomic positions x and momenta 
p after passing the interaction zone. For sufficiently large 
detuning between the field and atomic frequencies, one 
can adiabatically eliminate the excited-state amplitude of 
an atom B and derive a dimensionless effective Hamil- 
tonian for the center-of-mass motion 

o 

H e g= — e(l — sin (3t) sin x, (2) 

where x = kf < x > and p = < p > / hkf are the scaled 
expectation values for the atomic position and momen- 
tum, uj r — Hkf/mQo and /3 the scaled recoil and modu- 
lation frequencies, e and r = Hot the scaled well depth 
and time, respectively. 

Solving the respective Hamilton equations we compute 
outgoing atomic positions x ut at r = 500 as a function 
of ingoing positions x with uj r = 0.001, e = 0.625, and 
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FIG. 1: The outgoing positions of far-detuned atoms s out 
scattered by an amplitude-modulated standing wave as a 
function of initial positions xo in units kj 1 . 
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[3 = 0.07. Atoms in the beam are supposed to be ini- 
tially in the ground state, uniformly distributed along 
the re-axis with the same value of the initial momentum 
Po = 40.1 and strongly detuned with 5 = (ujf — uj a ) /^o = 
16. Chaotic atomic transversal motion in the modulated 
standing wave results in a self-similar structure of the 
scattering function which reveals itself under successive 
magnifications (see FIG. |l|). Let us consider the initial in- 
terval of xq (that may be identified with a cavity length) 
as a scattering region. Chaos implies that there exist 
atomic trajectories that never escape this region, and so 
there exist atoms never leaving the cavity. A set of ini- 
tial positions corresponding to nonescaping trajectories 
forms a fractal. We have found singularities computing 
the time atoms need to go out off the initial interval (to 
leave the cavity) . Similar results with appearing fractals 
have been obtained with far-detuned atoms scattered by 
a phase-modulated standing wave. Fractal-like peculiari- 
ties should appear in standard experiments on scattering 
atoms at standing waves fiofl in the case of their ampli- 
tude or phase modulations. 

3. Fractals in the atom-field dynamics. We wish 
to illustrate a generic effect of fractals in cavity QED and 
use the full Hamiltonian ([j]) without any modulation, i.e. 
f(t) = 1. A comparatively large average number of pho- 
tons is supposed in a single-mode cavity to sustain atom- 
field Rabi oscillations. By using the following expecta- 



tion values: 



x, p, 



= < er 2 >, u 



< a'(T_ + a<T+ >, and 



v = i < (vb- — air + >, we translate the Heisenberg equa- 
tions with the Hamiltonian ([!]) into the closed nonlinear 
semiclassical system 11 



x = u> r p, p = — usinx, u = Sv, z = — 2i>cosa 
v = -Su + [(2N - \)z - 3z 2 /2 + 1/2] cosx, 



(3) 



where dot denotes a derivative with respect to r. The di- 
mensionless control parameters are the recoil frequency 
uj r , the detuning 5, and the number of excitations 
N =< a) a + (o z + l)/2 >. The integral of motion, 



FIG. 2: (a) The schematic diagram showing a standing-wave 
microcavity with detectors, (b) Sample atomic trajectories on 
the plane of atomic momentum p (in units hkf) and position 
x (in units of kj 1 ). 



E = uj r p 2 /2 — ucosx — Sz/2, reflects the conservation 
of energy. At exact resonance, 5 = 0, the energy of the 
atom-photon interaction u is conserved, and the system 
(^) becomes integrable with solutions describing regular 
atomic center-of-mass motion in a potential well or over 
potential hills and atom-field Rabi oscillations modulated 
by the translational motion Out of resonance, the 
atomic translational motion is described by the equation 
for nonlinear parametric pendulum, 'x + lo t u(t) sin a; = 0, 
that has been analytically shown |ll]] to produce weak 
chaos even in the case of the simplest harmonic modula- 
tion u{t) ~ cos fir caused by the Rabi oscillations with a 
constant frequency Q = VS 2 + 4iV. In fact, the Rabi os- 
cillations is an amplitude- and frequency-modulated sig- 
nal which may induce pronounced erratic motion of the 
atomic center-of-mass inside the cavity. The motion is 
very complicated nearby the unperturped separatrix of 
the nonlinear pendulum where the period of oscillations 
goes to infinity. In this zone small changes in frequency, 
caused by respective small changes in energy, may lead 
to dramatic changes in phase which are the ultimate rea- 
son of exponential instability of atomic motion in a pe- 
riodic standing wave. The speculations above have been 
confirmed in our numerical experiments ^ |9| where posi- 
tive values of the maximal Lyapunov exponent have been 
found in the following ranges of the control parameters: 
u r > 10~ 3 , N < 10 2 , and | 5 | < 1.5. These magnitudes 
seem to be reasonable with available optical microcavi- 
ties in which the strength of the atom-field coupling may 
reach the values of the order of flo ~ 2-7T • 10 8 Hz, and 
the conditions of strong coupling are fulfilled for both the 
internal and external degrees of freedom ^ [| . 

Fractals in the dependencies x ou t(a;o) andp ou t(po) with 
atoms being transversely injected into a stationary stand- 
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FIG. 3: Atomic fractal in the strongly coupled atom-field sys- 
tem with three degrees of freedom. Exit time of atoms T in 
units O.Q vs the initial atomic momentum po in units hkf. 



ing wave have been found with Eqs. (||). In FIG. |^a we 
depict the scheme of another gedanken experiment with a 
Fabry-Perot optical microcavity and two detectors that 
may lead to another type of atomic fractals in cavity 
QED. To avoid complications that are not essential to 
the main theme of this work, we consider the cavity with 
only two standing-wave lengths. Atoms, one by one, are 
placed at the point x — with different initial values 
of the momentum poi along the cavity axis. We mea- 
sure the time when an atom reaches one of the detectors, 
the exit time T, and study the dependence T(p ) under 
the other equal initial conditions imposed on the atom 
and the cavity field. At exact resonance, optical poten- 
tial coincides with the standing-wave structure, and the 
analytical expression for the dependence in question can 
be easily found to be the following: T(po) = °o if po < 
Pcr/V2, P = Per; T(p ) = F[a,icsm(-l/K\/2),K]/y/u r uo 
if Pcr/\/2 ^ po < p CI :, and T(p ) = 2 F[37r/4, K]/uj r po 
if po > p CI . Here K = po^uj r /uo/2 is the modu- 
lus of the first-order incomplete elliptic integral F and 
p cr = 2^/ uo/uJr is the amplitude value of the atomic mo- 
mentum on the separatrix. Atoms with po < p CI /V2 are 
trapped in a potential well. Atoms with po > p cr fly 
through the cavity in one direction and are registered by 
one of the detectors. 

Out of resonance (S ^ 0), the atomic motion has 
been numerically found ^) and analytically proven 
Jrl"| to be chaotic with positive values of the maximal 
Lyapunov exponent. FIG. || shows the function T(p ) 
with the normalized detuning 5 = 0.4, the recoil fre- 
quency uj r = 1CP 3 , the number of excitations N = 10, 
and the following initial conditions: xq = 0, zo = 0, 
and wo = i>o = 2.17945. Atoms with comparatively 
small initial momentum, pq < 63, cannot reach the de- 
tectors because they are trapped in a well. With in- 
creasing po, the exit time function demonstrates an in- 
termittency of smooth curves and complicated structures 
that cannot be resolved in principle, no matter how large 



the magnification factor. FIG. [p shows magnification 
of the function for the small interval 86.7 ^ po ^ 88.3. 
Further magnifications reveals a self-similar fractal-like 
structure that is typical for Hamiltonian systems with 
chaotic scattering jl^, The length of the T(po) func- 
tion L = | Ti + \ — Ti | depends on the interval e 
of the partitioning of the momentum interval po- Wc 
compute the length of the fractal curve in FIG. || to 
be L(e) ~ e" 1 , where 7 ~ 0.84 ± 0.02 is a fractal di- 
mension simply connected with the Hausdorff dimension 
d= l+7~ 1.84 ±0.02. 

The exit time T, corresponding to both smooth and un- 
resolved po intervals, increases with increasing the mag- 
nification factor. It follows that there exist atoms never 
reaching the detectors in spite of the fact that they have 
no obvious energy restrictions to leave the cavity. Tiny 
interplay between chaotic external and internal dynam- 
ics prevents these atoms from leaving the cavity. The 
similar phenomenon in chaotic scattering is known as dy- 
namical trapping. Different kinds of atomic trajectories, 
which are computed with the system (J3J) , are shown in 
FIG. |2|b. A trajectory with the number to transverses 
the central node of the standing-wave before being de- 
tected m times and is called m-th trajectory. There 
are also special separatrix-like mS-trajectories following 
which atoms in infinite time reach the stationary points 
x = ±7m (n = 0,1,2,...), p — 0, transversing m times 
the central node. These points are the anti-nodes of the 
standing wave where the force acting on atoms is zero. In 
difference from separatrix motion in the resonant system 
(5 = 0) with the initial atomic momentum p cr , a de- 
tuned atom can asymptotically reach one of the station- 
ary points after transversing the central node m times. 
The smooth po intervals in the first-order structure in 
FIG. |^a correspond to atoms transversing once the cen- 
tral node and reaching the right detector. The unresolved 
singular points in the first-order structure with T = 00 at 
the border between the smooth and unresolved po inter- 
vals are generated by the 1 ^-trajectories. Analogously, 
the smooth and unresolved po intervals in the second- 
order structure in FIG. [|b correspond to the 2-nd or- 
der and the other trajectories, respectively, with singular 
points between them corresponding to the 2S'-trajectories 
and so on. 

There are two different mechanisms of generation of 
infinite detection times, namely, dynamical trapping 
with infinite oscillations (to = 00) in a cavity and the 
separatrix-like motion [m 7^ 00). The set of all initial 
momenta generating the separatrix-like trajectories is a 
countable fractal. Each point in the set can be specified 
as a vector in a Hilbert space with to integer nonzero com- 
ponents. One is able to prescribe to any unresolved inter- 
val of m-th order structure a set with to integers, where 
the first integer is a number of a second-order structure 
to which trajectory under consideration belongs in the 
first-order structure, the second integer is a number of a 
third-order structure in the second-order structure men- 
tioned above, and so on. Such a number set is analogous 
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FIG. 4: Exit time distributions in the regimes of strong chaos 
with S = 0.4 (crosses) and weak chaos with 5 — 0.05 (cir- 
cles) with the same initial conditions and the other control 
parameters, as in FIG. H. 



to a directory tree address: " <a subdirectory of the root 
directory>/<a subdirectory of the 2-nd level>/<a sub- 
directory of the 3-rd level>/...". Unlike the separatrix 
fractal, the set of all initial atomic momenta leading to 
dynamically trapped atoms with m = oo seems to be 
uncountable. 

The dependence of the maximal Lyapunov exponent 
of the set (|^) on the detuning S have been found almost 
the same as compared with a slightly different version 
of (H) considered in ||. It shows a deep gorge around 
5 = 0, maxima around | 5 | ~ 0.5, and falls to zero at 
| 5 | > 1.5. We collect an exit time statistics, by count- 
ing atoms with different initial momenta reaching the 
detectors, in the regimes of comparatively strong chaos 
(5 = 0.4) and weak chaos (5 = 0.05). The plots of the 
respective histograms of exit times are shown in FIG. ||. 
The probability distribution function P(T) with almost 
resonant atoms decays rapidly and demonstrates a single 



local maximum around T ~ 140. In the regime of strong 
chaos, P(T) demonstrates a few local peaks (the first one 
around T ~ 140) and a long tail up to T ~ 500. 

4. In summary, fractals in single-atom standing-wave 
cavity QED typically arise in the nonlinear models of the 
atom-field interaction which are chaotic in the semiclas- 
sical limit. It has been shown as for the non-autonomous 
Hamiltonian system (||) with one and half degree of free- 
dom describing the nonlinear motion of a "structure- 
less" atom in a far detuned amplitude (phase)-modulated 
standing-wave field, as for the autonomous system (||) 
with three degrees of freedom including internal and ex- 
ternal atomic and field ones. The fractals of the types to 
be found in this Letter should appear also in the center- 
of-mass motion of an ion confined in a Paul trap and in- 
teracting with a laser field Il4l] . Different effective Hamil- 
tonians have been introduced to describe this interaction. 
The simplest one [O] is of a periodically driven linear os- 
cillator which is chaotic in the classical limit |l6| . 

We have considered the semiclassical chaotic scatter- 
ing. In quantum dynamics there are no individual atomic 
trajectories. It implies lacking of singularities in scatter- 
ing functions, but classical fractals should manifest them- 
self in fluctuations of corresponding quantum S-matrices 
whose structures for generic quantum scattering models 
have been shown to be described by random matrices 
fl7f . We plan to treat the quantum-classical correspon- 
dence in chaotic scattering of atoms by a standing wave 
in the future. 
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